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Abstract: Information on rare earth elements (REEs) in soils and plants of the Qinghai-Tibet Plateau is 
very limited. Therefore, in this study, we performed field sampling to explore the geochemical signatures 
and human health risk of REEs in soils and plants of the northeastern Qinghai-Tibet Plateau, China. A 
total of 127 soil samples and 127 plant samples were collected from the northeastern Qinghai-Tibet 
Plateau to acquire the geochemical signatures and related human health risks of REEs. The mean total 
concentrations of REEs in soils and plants of the study area reached 178.55 and 10.06 mg/kg, 
respectively. The light REEs in soils and plants accounted for 76% and 77% of the total REEs, 
respectively. REEs showed significantly homogenous distribution in soils but inhomogeneous distribution 
in plants of the study area. Characteristic parameters indicated that light REEs were enriched and 
fractionated significantly, while heavy REEs were moderately fractionated in soils and plants. REEs in 
soils and plants showed significantly negative Europium anomaly. Cerium showed slightly positive 
anomaly in plants and slight anomaly in soils. The normalized distribution patterns of REEs were 
generally similar in the analyzed soils and the corresponding plants of the study area. The average 
bio-concentration factor of REEs ranged from 0.0478 (Scandium) to 0.0604 (Europium), confirming a 
small accumulation of REEs by plants. Health risks caused by REEs in soils and plants were negligible, 
while risks for adults were lower than those for children. This study provides important information on 
КЕЕ in soils and plants of the northeastern Qinghai-Tibet Plateau. 


Keywords: rare earth elements; geochemical signatures; human health risk; carcinogenic risk; bio-concentration factor; 
Qinghai-Tibet Plateau 


Citation: LI Leiming, WU Jun, LU Jian, ZHANG Xiying, XU Juan. 2022. Geochemical signatures and human health risk 
evaluation of rare earth elements in soils and plants of the northeastern Qinghai-Tibet Plateau, China. Journal of Arid Land, 
14(11): 1258-1273. https://doi.org/10.1007/s40333-022-0107-8 


1 Introduction 


Rare earth elements (REEs) encompass the yttrium (Y), the elements of the lanthanide series 
(atomic number from 57 to 71), and scandium (Sc), which exhibit similar environmental 


"Corresponding author: WU Jun (E-mail: wujunlisa@ 163.com) 
Received 2022-07-06; revised 2022-10-09; accepted 2022-10-14 


© Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Science Press and Springer-Verlag GmbH Germany, 
part of Springer Nature 2022 


http://jal.xjegi.com; www.springer.com/40333 


LI Leiming et al.: Geochemical signatures and human health risk evaluation of rare... 


behaviors (Ferreira et al., 2021). Promethium (Pm) does not exist in the Earth's crust since the 
radioactive decay of Pm is fast (Migaszewski and Gatuszka, 2015). REEs generally exist together 
in natural environment due to similar physical and chemical properties (Groenenberg et al., 2010; 
Khan et al., 2017). Many studies reported that pedogenic processes, soil characteristics, and 
parent material affected REEs background contents in soils (Yao et al., 2010; Wu et al., 2019). 
КЕЕ; are not rare from a geochemical perspective because the abundances of REEs in soils are 
similar to those of zinc (Zn), lead (Pb), and copper (Cu) and higher than those of mercury (Hg), 
Stannum (Sn), cobalt (Co), and silver (Ag) (Wang et al., 1998). 

Heavy REEs (ZHREEs; the sum Y, Sc, and from gadolinium (Gd) to lutetium (Lu)) and light 
REEs (ZLREEs; the sum from lanthanum (La) to europium (Eu)) compose REEs based on the 
mass and atomic numbers of REEs. The concentrations of XHREESs are generally lower than the 
concentrations of XLREEs (Henderson, 1984; Durn et al., 2021). The XLREEs are more mobile 
and available for plant absorption than XHREEs (Brioschi et al., 2013). The concentrations of 
XLREEs in leachate are higher than those of XHREEs so that ће XLREEs are more available for 
root adsorption than XHREEs (Aide and Aide, 2012). The bioavailability of REEs in soil is a 
considerable factor for the accumulation of REEs by plants (Li et al., 2013; Liu et al., 2018). The 
concept of bioavailability also was defined the amount of human adsorption in recent years (Li et 
al., 2018; Ferreira et al., 2022; Wang et al., 2022). 

The global annual demands of REEs have increased significantly from 3.00x10* to 1.35x10? t 
during 1980—2015 (Barta, 2007; Huang et al., 2016). China is the largest producer of REEs in the 
world, accounting for approximately 9796 of the global supply. The largest REEs ore mine is 
located in the northern China (Tang et al., 2020). REEs have been widely used in medical 
technology, industry, and agriculture (Hu et al., 2006; Pagano et al., 2015; Wang and Liang, 2016; 
Ou et al., 2022). Various field activities have released lots of REEs into the environment (Zhang 
et al., 2014; Pagano, 2016). Distribution, environmental behaviors, and removal of different 
contaminants, such as trace elements (Zheng et al. 2010a, b; Liang et al, 2013; 
Rodríguez-Barranco et al., 2014; Lu et al., 2021; Pirarath et al., 2021; Cui et al., 2022; He et al., 
2022), antibiotics and corresponding resistance genes, perfluorooctane sulfonate, polycyclic 
aromatic hydrocarbons, and microplastics, have attracted attention widely (Liu et al., 2019; Fu et 
al., 2021; Qiao et al., 2021; Lu et al., 2022a, b). However, studies regarding the negative effect of 
REEs on natural ecosystems and humans are limited. REEs accumulate in various environmental 
media (including soil, atmosphere, and water body) after entering into the environment, and 
finally influence ecosystem via human activities or human health via REEs accumulation in the 
human bodies (Mihajlovic et al., 2015; Pagano, 2016; Wang et al., 2017; Wang et al., 2019; 
Godwyn-Paulson et al., 2022). The behavior of REEs in natural environment is strongly affected 
by soil redox conditions, organic matter content, pH, weathering processes, climate conditions, 
and parent rock (Zhang et al., 2014; Pagano, 2016; Allajbeu et al., 2016; Huang et al., 2019). 
There are three ways including ingestion, dermal adsorption, and inhalation for entering the 
human body (Meryem et al., 2016; Guo et al., 2019). Previous studies have proven that REEs 
accumulated in the human body can cause numerous diseases (Zhang et al., 2000; Zhao et al., 
2017; Ferreira et al., 2022). 

Soil is an important medium for receiving pollutants and nutrients to play an important role in 
ecosystems. In recent decades, the geochemical characteristics and health hazards of REEs have 
attracted much attention (Kan et al., 2017; Huang et al., 2019; Guo et al., 2019; Wu et al., 2019; 
Malhotra et al., 2020). However, information about distribution patterns, accumulation in soils 
and plants, and health risk evaluation of REEs is quite scarce in the study area. The main 
objectives of this study were to (1) investigate the concentration and distribution patterns of REEs 
in soils and plants of the northeastern Qinghai-Tibet Plateau; (2) calculate bioavailability and 
mobility characteristics of REEs; and (3) evaluate the potential human health risk of REEs. 
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2 Materials and methods 


2.1 Study area and soil and plant samples 


The study area is situated in the northeastern Qinghai-Tibet Plateau of China, which spans 
4.13x10° km? (90°37'-101°05'E, 35?01'-38?43'N). The study area has an average annual 
precipitation of 50-450 mm and an annual average temperature of —5.7—8.5?C. A total of 127 soil 
samples (soil depth: 0—20 cm) and 127 plant samples were collected from June 14 to June 29 of 
2017 by following standard sampling method (State Environmental Protection Administration, 
2004) (Fig. 1). All samples were placed in plastic bags, transported back to the laboratory, and 
stored at -80?C. Soil samples were prepared by following the study conducted by Li et al. (2018). 
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Fig.1 Study area and soil and plant sampling sites 


2.2 Analysis of rare earth elements (REEs) 


The soil samples were air-dried, crushed, and then sieved using 0.074-mm sieves. The soil 
samples were calcined at 700?C for 2 h to eliminate the organic matter. According to related 
standard (State Environmental Protection Administration, 2004), 100 mg of each soil power was 
added along with 3 mL HNO; with high purity, 2 mL НСІ with high purity, and 2 mL HF with 
high purity in a 100-ml Teflon digestion tube (Wu et al., 2019). These tubes were placed in a 
microwave dissolving system for 1 h and cooled. The cooled solution was dried at 180°C with 2 
mL of high quality HClO, and 15 mL of high purity HNO3. 

The plant (Potentilla anserina L.) samples were pre-treated by following the 
previously-reported procedures (Shen et al., 2018). All parts of plant samples were washed 
thoroughly to remove soil particles, washed five times with ultra-pure water, blotted extra water 
with tissue paper, and dried. Then the samples were grinded to fine powder for digestion. 
Approximately 100 mg of dry plant powder was digested with 3 mL НСІ (37%), 5 mL HNO3 
(65%), and 1 mL HF (65%) (Ayrault et al., 2001; Li et al., 2018). The other digestion procedures 
for plant samples were same as soils. Using testing instrument (Agilent7900, Agilent Inc, Santa 
Clara, USA), we measured the concentration levels of La, cerium (Ce), praseodymium (Pr), 
neodymium (Nd), samarium (Sm), Eu, Gd, terbium (Tb), dysprosium (Dy), holmium (Ho), 
erbium (Er), thulium (Tm) ytterbium (Yb), Lu, Y and Sc in samples. 


2.3 Distribution characteristics of REEs 


North American shale composite (NASC) was used to normalize REEs concentrations to 
eliminate the "Oddo-Harkin effect" of REEs (Taylor and Mclennan, 1985). Normalization values 
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were calculated by the ratio of the REEs concentrations to the corresponding element 
concentrations of NASC in this study. The parity effects of REEs were eliminated after the 
normalized REEs by NASC (Chen and Yang, 2010; Wu et al., 2019). 

Characteristic parameters which indicated the distribution of REEs in soils and plants included 
fractionation between XLREEs and XHREEs ((La/Yb)n), fractionation between XLREEs 
((La/Sm)y), fractionation between XHREEs ((Gd/Yb)n), Ce anomaly (dce), and Eu anomaly (ёва); 
we calculated these parameters based on normalized values of NASC (Taylor and 
Mclennan,1985; Mclennan, 1989; Chen and Yang, 2010; Wu et al., 2019; Yu et al., 2019), the 
calculation process is as follows: 


Га, „у. Ла 
(La/Sm),, = ( sample asc) : (1) 
(ЗИ гааг /SMyasc) 
(Gd/Yb) y= (Gd sample Gd wAsc) (2) 
(ҰБ, / Үд хдас) | 
(La/Yb) N= (La sample /La NASC ) (3) 
СҮН, ole [Yb NASC ) 
M (Се ipis /Cexasc) (4) 
Al (La sample Ла nasc) x (Pr, Т.атрїе /Ptyasc) 
5. = (La sample’ La nase) (5) 
Eu~ , 
(STA, cols /SMyasc)/ (Gd sample /Gdyasc) 
where Lasample (mg/kg), паре (mg/kg), Gdsample (mg/kg), YDsample (mg/kg), Cesample (mg/kg), and 


Prsampie (mg/kg) are the contents of La, Sm, Gd, Yb, Ce, and Pr in soil or plant samples, respectively; 
and Lawasc (mg/kg), Smyasc (mg/kg), Gdyasc (mg/kg), Ybyasc (mg/kg), Cenasc (mg/kg), and 
Pryasc (mg/kg) are the contents of La, Sm, Gd, Yb, Ce, and Pr in reference material, respectively. 
The content of REEs in reference materials can be found in Taylor and Mclennan (1985). 

2.4 Bio-concentration factor (ВСЕ) 

Bio-concentration factor (BCF) is usually adopted to calculate the potential transfer of various 


contaminants from soils to plants (Liu et al., 2015; Ingrid et al., 2022; Nicolas et al., 2022). This 
approach was also employed to examine REEs. The BCF is calculated as follows: 


X 
BCF =, (6) 
soil 
where BCF; is the BCF value of the i target REEs; and Xpian (mg/kg) and X;oii (mg/kg) are the 
concentrations of the target REEs in soil and plant samples, respectively. 


2.5 Human health risk evaluation 


Three exposure pathways including ingestion, inhalation, and dermal adsorption of REEs in soils 
and plants might exert threat to human health (Zheng et al., 2020). This study used formulas from 
the Environmental Protection Agency of USA and other articles (USEPA, 1996; Chen et al., 
2022; Li et al. 2022). The hazard quotients posed through three ways can be calculated using the 
following equations: 


НО, = х IngR x EFx ED 1056, (7) 
е "BWxATxCFD 
НО, - Cx InhR x EFx ED (8) 


PEF x BW x AT x CFD 
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HG oq AA TAL qu. (9) 
BW x AT x CFD 
HI=HQing + HQ; + HQ... (10) 


where HQing, HQinn, and HQae are the hazard quotient of ingestion, inhalation, and dermal 
adsorption, respectively; C is element concentration (mg/kg); BW is body weight (kg); EF is 
exposure frequency (d/a); ED is exposure duration (a); AT is averaging time (d); CFD is 
corresponding reference dose (mg/(kg-d)); IngR (mg/d) and InhR (m?/d) are ingestion rate and 
inhalation rate, respectively; PEF is particle emission factor (m?/mg); SAF is skin adherence 
factor (mg/(cm?.d)); SAE is skin area exposed (cm?); DAF is dermal absorption factor; and HI 18 
the sum of HQing, HQinh, and HQaer. 

Reference values used in the equations are from previous studies (USEPA, 2002; Duan, 2012; 
Faiz et al., 2012; Chen et al., 2022; Li et al., 2022). We used a unique reference dose value (0.02 
mg/(kg-d)) for all REEs in this study (Li et al., 2013; Sun et al., 2017). The values of hazard 
quotient above 1.0 are considered non-carcinogenic risk to human health and the values of hazard 
quotient below 1.0 are considered safe for the human body (Wu et al., 2018; Li et al., 2022). 
Carcinogenic risk was calculated by lifetime mean daily dose (LMDD, mg/(kg-d)) (Sun et al., 
2017; USEPA, 2004) and used the following formulas: 


LMDD- CxEF " | CR баа X ED enia iE CR аа X ED aut | aD 
AT BW anita BW dutt 
Carcinogenic risk = LMDD*xSF, (12) 


where CReniia and CRadut are contact rate for child and adult, respectively; EDcnia and Ер аш are 
exposure duration for children and adults (a), respectively; BW.nia and BWaau are body weight of 
children and adults (kg), respectively; and SF is slope factor, which is 3.2х10:12 in the study (Sun 
et al., 2017). The unacceptable risk for regulatory purposes is higher than 1.0x10~ (Faiz et al., 
2012; Chen et al., 2022). 


2.6 Statistical analysis 


The statistical data of REEs were obtained by Excel while the correlation of REEs with soils or 
plants was studied by Pearson correlation analysis by using SPSS 22.0. Digital terrain model was 
established by applying inverse distance weight interpolation (Li et al., 2013; Xu et al., 2016; 
Wang et al., 2019; Wu et al., 2021). ArcGIS 10.2 processed corresponding distribution for data on 
human health risk values. 


3 Results and discussion 


31 КЕЕ; of soils 


Detailed statistical results on the concentrations of total REEs and individual REEs in 127 soil 
samples were showed in Table 1. The concentrations of REEs normalized NASC, background 
values in Qinghai Province of China, and background values in China were presented (Taylor and 
Mclennan, 1985; CNEMC, 1989; MEPC, 1990). The concentrations of REEs in soils ranged from 
55.29 to 306.80 mg/kg with an average of 178.55 mg/kg (Table 1). The concentration of REEs in 
soils from most of sampling sites were generally similar in the study area. The REEs exhibited 
significantly good homogeneity of distribution in soils. The total concentrations of REEs in 127 
soil samples exceeded 100.00 mg/kg. Two soil samples with the highest concentrations of REEs 
were located in an industrial area (more than 300.00 mg/kg). The results indicated that 
anthropogenic activities, particularly mining or industrial activities, may contribute slightly to the 
accumulation of REEs in some sites of study area. 

The average concentrations of REEs of NASC, background values in Qinghai Province of China, 
and background values in China were lower than the concentrations in soils of the study area 
(Taylor and Mclennan, 1985; CNEMC, 1990; MEPC, 1990). Likewise, the average concentrations 
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of XLREESs in soils (135.26 mg/kg) were lower than the background values in China (147.92 mg/kg) 
and higher than that of NASC (73.50 mg/kg) and background values in Qinghai Province of China 
(126.37 mg/kg). Mean concentrations of SHREEs (43.29 mg/kg) were higher than the XHREEs of 
NASC (33.40 mg/kg), background values in Qinghai Province of China (35.46 mg/kg), and 
background values in China (38.84 mg/kg) (Tablel). XLREEs in soils accounted for 76% of the 
total REEs, which was consistent with Tyler's results (Tyler, 2004). The results indicated that 
concentrations of XLREEs were higher than those of XHREES, similar to the results of NASC. The 
correlation analysis illustrated that there was an obvious positive correlation among REEs of soils 
(Fig. 51), while the relationship among Eu, Er, Sc, and other REEs was poor. These results indicated 
that most of REEs might have similar sources in the soils. 


Table 1 Statistic parameters of rare earth elements (REEs) in soils and plants of the northeastern Qinghai-Tibet 
Plateau, China 


Soils (mg/kg) Plants (mg/kg) Background 


Rare earth values in Background 
element NASC’ Qinghai o 
and Min Mean Max SD Min Mean Мах sp (mg/kg) Province of (soils) 
parameter China (soils) (mg/kg) 
(mg/kg) 
La 8.36 31.04 54.26 6.6 0.10 1.76 8.66 1.7 16.00 32.80 39.70 
Ce 17.09 63.66 112.43 13.2 0.44 3.78 16.70 3.4 33.00 58.30 68.40 
Pr 1.82 7.32 12.50 1.6 0.02 0.39 1.93 0.4 3.90 5.87 TAT 
Nd 6.95 26.95 46.38 5.8 0.08 1.43 7.19 14 16.0 23.70 26.40 
Sm 1.33 5.18 8.94 1.1 0.02 0.28 1.45 0.3 3.50 4.77 5.22 
Eu 0.33 1.11 3.26 0.3 0.00 0.07 0.32 0.1 1.10 0.93 1.03 
Gd 1.31 5.10 8.81 1.0 0.02 0.28 1.41 0.3 3.30 4.15 4.60 
Tb 0.18 0.68 1.24 0.1 0.00 0.04 0.19 0.0 0.60 0.60 0.63 
Dy 0.99 3.56 6.68 0.7 0.01 0.20 1.02 0.2 3.70 3.80 4.13 
Ho 0.20 0.69 1.29 0.1 0.00 0.04 0.20 0.0 0.80 0.77 0.87 
Er 0.64 2.19 6.17 0.7 0.01 0.11 0.58 0.1 2.20 2.32 2.54 
Tm 0.08 0.28 0.51 0.1 0.00 0.02 0.08 0.0 0.30 0.33 0.37 
Yb 0.53 1.82 3.32 0.4 0.01 0.11 0.52 0.1 2.20 2.07 2.44 
Lu 0.08 0.27 0.47 0.1 0.00 0.02 0.08 0.0 0.30 0.32 0.36 
Y 5.13 18.20 34.64 3.9 0.07 1.07 5.43 1.0 20.00 21.10 22.90 
Sc 4.86 10.51 21.93 32 0.03 0.47 2.48 0.4 30.00 10.42 / 
REEs 55.29 178.55 306.80 35.8 1.03 10.06 48.25 9.3 136.90 172.25 186.76 
XLREEs 35.88 13526 23442 282 0.84 7.71 36.27 7.2 73.50 126.37 147.92 
XHREEs 18.34 43.29 72.39 8.8 0.18 2.35 11.98 22 33.40 35.50 38.84 
MOREAT (275 414 — 600 04 295 420 674 06 220 3.56 3.81 
(La/Yb)n 1.71 2.36 3.61 0.3 1.52 2.22 3.01 0.3 / 2.18 2.24 
(La/Sm)n 1.16 1.31 1.54 0.1 1.08 1.36 1.67 0.1 / 1.50 1.66 
(Gd/Yb)n 1.50 1.87 2.51 0.2 1.29 1.78 3.06 0.2 / 1.34 1.26 
бсс 0.95 1.01 131 0.1 0.93 1.23 3.06 0.4 / 1.01 0.97 
Or. 0.49 0.67 2.02 0.1 0.51 0.81 3.14 0.3 / 0.65 0.65 


Note: La, lanthanum; Ce, cerium; Pr, praseodymium; Nd, neodymium; Sm, samarium; Eu, europium; Gd, gadolinium; Tb, terbium; Dy, 
dysprosium; Ho, holmium; Er, erbium; Tm, thulium; Yb, ytterbium; Lu, lutetium; Y, yttrium; Sc, scandium; REES, rare earth elements; 
XLREEs, light REEs (the sum from La to Eu); XHREEs, Heavy REEs (the sum Y, Sc, and from Gd to Lu); (La/Yb)ny, fractionation 
between XLREEs and XHREEs; (La/Sm)w, fractionation between XLREEs; (Gd/Yb)y, fractionation between XHREEs; ӧс, Ce anomaly; 
бв, Eu anomaly; Min, minimum; Max, maximum; SD, standard deviation; NASC, North American shale composite. * refers to Taylor 
and Mclennan (1985), ° refers to МЕРС (1990), * refers to CNEMC (1989), and / refers to no data. 
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3.2 REEs of plants and BCF 


Table 1 also showed the concentrations of total and individual REEs in plants. The order of 
concentrations of REEs in plants was basically the same as that in soils. The distribution of REEs 
in plants was extremely uneven in the study area. The average concentrations of REEs in plants 
were variable in different sites. The average concentration of total REEs was 10.06 mg/kg, and 
while the concentrations of total REEs ranged from 1.03 to 48.25 mg/kg in plants. The average 
concentration of XLREEs was 7.71 mg/kg in plants, and the concentration range of XLREEs was 
0.84 to 36.27 mg/kg. The concentrations of SHREEs ranged from 0.18 to 11.98 mg/kg, with an 
average value of 2.35 mg/kg. XLREEs covered 77% of total REEs in plants of the study area. La, 
Ce, Nd, and Y in plants accounted for more than 8046 of total REEs. The characteristics of REEs 
in soils of the study area were well correlated with the high fraction of the above REEs and the 
distribution order of REEs concentrations. The correlation analysis showed that there was a 
significantly positive correlation among the REEs in plants to prove that all REEs in plants came 
from the same sources (Fig. S1b). 

BCF can often reflect a plant's capacity to accumulate certain components. The statistical 
information on BCF values of REEs at various sites was shown in Figure 2. At the sampling sites, 
BCF values ranged from 0.004 to 0.289 for all target REEs. The average BCF values for REEs 
followed the order of Eu>Lu>Y>Ce>Yb>Ho>Tm>Dy>La>Tb>Er>Gd>Sm>Nd>Pr>Sc. The 
average BCF value of Eu in the study area was high compared to other REEs, indicating a 
relatively high mobility of transplants in soils compared to other REEs. In the same study 
location, the mean BCF values of REEs were significantly lower than those of heavy metals (Li et 
al., 2018). However, the average BCF of REEs was lower than that of heavy metals, indicating 
that REEs were more difficult to transfer from soils to plants (Hu et al., 2013; Jeelani et al., 2017; 
Liu et al., 2017; Li et al., 2018). BCF values of REEs in all sites were less than 1.000, indicating 
relatively low bioaccumulation. The results were similar to the recent studies (Ingrid et al., 2022; 
Nicolas et al., 2022; Tao et al., 2022). 
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Fig. 2 Bio-concentration factor (BCF) values of rare earth elements (REEs) in the study area. La, lanthanum; 
Ce, cerium; Pr, praseodymium; Nd, neodymium; Sm, samarium; Eu, europium; Gd, gadolinium; Tb, terbium; Dy, 
dysprosium; Y, yttrium; Но, holmium; Er, erbium; Tm, thulium; Yb, ytterbium; Lu, lutetium; Sc, scandium. The 
boxes represent the range from the lower quantile (Q25) to the upper quantile (Q75). The dots and horizontal lines 
inside the boxes represent the means and medians, respectively. The dots outside the boxes represent outliers. The 
upper and lower whiskers show the range within 1.5IQR (interquartile range). 
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3.3 REEs Characteristic parameter in soils and plants 


Table 1 and Figure 3 showed the characteristic parameters of REEs in soils and plants, including 
бии, ӧсе, XLREES/EZHREEs, (La/Sm)w, (Gd/Yb)n, and (La/Yb)w. The average value of 
XLREEs/ZHREESs in soils was 4.14, and the range of XLREEs/XHREESs in soils was 2.75—6.09. 
The results were similar to previous report (Wu et al., 2018). The range of SLREEs/XHERRs in 
plants was 2.95-6.74 with an average value of 4.21, indicating that the values of XLREEs were 
more abundant in plants than soils. The characteristic parameter of (La/Sm)w reflected the 
fractionation degree of XLREEs. The average value of (La/Sm)n in soils was 1.31, while the 
range of (La/Sm)n values was 1.16-1.54 in soils. The average value of (La/Sm)n in plants was 
1.36, while the range of (La/Sm)y values was 1.08-1.67. These results indicated that XLREEs in 
soils and plants were significantly fractionated in the study area. 
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Fig. З Concentrations of light REEs (ZLREEs) and heavy REEs (ZHREEs) as well as the ratio of XLREEs to 
XHREEs (ZLREEs/ZHREEs) in soil samples (a) and plant samples (b). Р-1-Р-127 represent the serial number of 
sampling sites. 


The characteristic parameter of (Gd/Yb)n illustrated the fractionation degree of XHREEs. The 
range of (Gd/Yb)n values in soils was 1.50—2.51 with an average value of 1.87, while (Gd/Yb)n 
values in plants ranged from 1.29 to 3.06 with an average value of 1.78 in the study area, 
illustrating that XHREEs in soils and plants had moderate fractionation. The characteristic 
parameter of (La/Yb)n showed the fractionation degree of XLREEs and XHREEs. The range of 
(La/Yb)n values in soils ranged from 1.71 to 3.61 with an average value of 2.36, while the range 
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of (La/Yb)w values in plants was 1.52-3.01 with an average value of 2.22, illustrating that 
XLREEs-enriched and the curve of the figure was right oblique (Fig. 4). These characteristic 
parameters in this study were similar to the background values in Qinghai Province of China and 
the background values in China. 
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Fig. 4 (a) normalization of average REEs values and background composition in soils; (b), normalization of 
average REEs values in plants. NASC, North American shale composite; BC, background composition. 
Soil/NASC represents normalization of average REEs values in soils, BC/NASC denotes normalization of 
background composition in soils, and Plant/NASC represents normalization of average REEs values in plants. 


Og, and doce generally indicate anomalies of Eu and Ce (Yao et al., 2010; Gill et al. 2018). The 
elements are positive anomalies when dgu or óc. values are greater than 1.05, conversely, they are 
negative anomalies when бєх or бсе values are lower than 0.95 (Zhao et al. 2017; Gill et al. 2018; 
Wu et al., 2019). The range of бв. values in soils was 0.49-2.02 with an average value of 0.67 in 
the study area. The average бв. value in plants was 0.81, ranging from 0.51 to 3.14. These results 
indicated Eu in soils and plants of the study area was significantly negative anomalies. The range 
of dce values in soils was 0.95-1.31 with an average value of 1.01. The range of óc. values in 
plants was 0.93—3.06 with an average value of 1.23. These results illustrated that Ce did not show 
significantly anomaly in soils while it showed slightly positive anomaly in plants. The gu and dce 
of REEs in soils and plants of the study area were very similar to the background values in 
Qinghai Province of China and the background values in China. 

Because REEs with even atomic numbers are always more numerous than those with odd 
atomic numbers, it is challenging to compare the abundance of REEs directly. The concentrations 
of REEs were normalized using NASC to remove "Oddo-Harkin effect" in the study area (Fig. 4). 
The normalization model can be adopted to compare the REEs abundance in various 
environmental media and to determine whether there is a shortfall or enrichment of REEs (Wu et 
al., 2019; Ferreira et al., 2021; Lin et al., 2021). The NASC-normalized distribution patterns of 
REEs were usually similar in soils and plants. The average REEs patterns of the soils and plants 
matched the average values patterns of the background values in Qinghai Province of China (Fig. 
4a and b). The distribution patterns of REEs in soils and plants of the study area were similar to 
the REEs patterns of the sediments and soils (Wu et al., 2019; Costa et al., 2021). Similar 
distribution patterns of REEs demonstrated that they primarily came from natural sources. 


3.4 Health risk of soils and plants 


Non-carcinogenic effects of REEs in soils and plants of the study area on human beings were 
calculated (Fig. 5; Tables S1 and S2). The total non-carcinogenic risks of REEs in soils for adults 
ranged from 0.0081 to 0.0450, while those in plants for adults ranged from 0.0002 to 0.0071. The 
range of total non-carcinogenic risks from REEs in soils for children was 0.0187 to 0.1038, while 


202211.00147v1 


chinaXiv 


ChinaXiv& f'ERBTII 


LI Leiming et al.: Geochemical signatures and human health risk evaluation of rare... 


the total non-carcinogenic risks of REEs in plants ranged from 0.0003 to 0.0163. The higher 
non-carcinogenic risks occurred in the eastern study area where soils and plants were underlain by 
human activities (Fig. 5). The lowest non-carcinogenic risk value was acquired from inhalation, 
followed by dermal adsorption and ingestion pathways for children and adults from soils and 
plants. Total non-carcinogenic risk values of soils and plants obtained in this study were less than 
1.0000, illustrating low non-carcinogenic risk of REEs. The non-carcinogenic risk posed by REEs 
from soil and plant ingestion were relatively low in the study area by considering the low REEs 
concentrations in soils and plants. 
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Fig. 5 Non-carcinogenic risk of REEs from soils and plants in the northeastern Qinghai-Tibet Plateau, China. (a 
and b), non-carcinogenic risk of REEs for adults and children from soils, respectively; (c and d), non-carcinogenic 
risk of REEs for adults and children from plants, respectively. HI is the sum of the hazard quotient of ingestion, 
inhalation, and dermal adsorption. 


Previous studies on the health evaluation of REEs in polluted areas have shown that the 
non-carcinogenic risk values exceeded more than 3.000x10^ for adults and 1.800х10:2 for 
children (Pagano et al, 2015; Ferreira et al., 2022), significantly higher than those values observed 
in plants from this study. The average non-carcinogenic risk values of La, Ce, and Nd in soils for 
adults was higher than 3.000x10?, while the non-carcinogenic risk value of Ce in soil for 
children was higher than 1.800х1072. Normally, children are more susceptible to ingesting REEs 
due to hand and mouth movements. In addition, children exhibited greater health risks for 
incomplete formation of physiologic system and low body weight (Sun et al., 2017; Li et al., 
2022). Based on the results of hazard quotient in Chinese cities, total non-carcinogenic risk values 
higher than 0.1000 for children should be more attention (Sun et al., 2017; Guo et al., 2019). 

The effect of carcinogenic risk of REEs in soils and plants of the study area on human beings 
were analyzed (Fig. 6; Tables S1 and S2). The total carcinogenic risk values of REEs in soils 
ranged from 1.593х107!9 to 8.839х107!9 (Fig. ба), while the total carcinogenic risk values of 
REEs in plants ranged from 2.976x10-? to 1.390х10:10 (Fig. 6b). Carcinogenic risk posed by 
REEs in the study area had similar distribution. Carcinogenic risk through ingestion, dermal 
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adsorption, and inhalation were all below safe level (<1.000x10~) in the study area, illustrating 
that carcinogenic risk generated through REEs to human body was acceptable in the study area. 
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Fig. 6 Carcinogenic risk (CR) of REEs from soils (a) and plants (b) in the northeastern Qinghai-Tibet Plateau, 
China 


4 Conclusions 


The distribution of REEs in soils was noticeably homogeneous although that in plants was not 
homogeneous. XLREEs in soils and plants account for 76% and 77% of total REEs, respectively. 
Most of REEs might have similar sources in soils or plants. The distribution patterns of REEs in 
soils and plants from the study area were generally similar. XLREEs had significant fractionation 
while XHREEs had moderate fractionation in soils and plants. XLREEs enriched in soils and 
plants. Significant negative anomaly of Eu occurred in soils and plants. Ce was not significantly 
anomaly in soils and had slightly positive anomaly in plants. The average BCF values followed 
the order of Eu>Lu>Y>Ce>Yb>Ho>Tm>Dy>La>Tb>Er>Gd>Sm>Nd>Pr>Sc. BCF values for 
REEs indicated that they had relatively low bioaccumulation. Eu had relatively high mobility of 
transplants of the study area. These characteristic parameters of REEs in soils and plants of the 
study area were similar to the background values in Qinghai Province of China and the 
background values in China. Both non-carcinogenic and carcinogenic risks posed by REEs in 
soils and plants were acceptable. Adults were less sensitive to REEs ingestion compared to 
children. 
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Fig. S1 Pearson correlations among rare earth elements in soils (a) and plants (b) of the northeastern 
Qinghai-Tibet Plateau, China. La, lanthanum; Ce, cerium; Pr, praseodymium; Nd, neodymium; Sm, samarium; 
Eu, europium; Gd, gadolinium; Tb, terbium; Dy, dysprosium; Ho, holmium; Er, erbium; Tm, thulium; Yb, 
ytterbium; Lu, lutetium; Y, yttrium; Sc, scandium. 


Table S1  Non-carcinogenic and carcinogenic risk for rare earth elements (REEs) in soils of the northeastern 
Qinghai-Tibet Plateau, China 


Non-carcinogenic risk for children Non-carcinogenic risk for adults Carcinogenic risks 
REE 
Ingestion Inhalation Dermal Total Ingestion Inhalation Dermal Total Ingestion Inhalation Dermal Total 
adsorption adsorption adsorption 


La  L02x10? 7.50х10% 2.86x10* 1.05x10? 4.37x10? 3.22х10% 1.75x10* 4.55x10? 8.00х107" 6.93x10 ? 2.52x10 ? 8.94х107! 
Ce  2.00x10? 1.54х10° 5.86х10* 2.15x10? 8.97x10? 6.60x10° 3.58х10* 9.33x10? 1.64х1079 1.42x10"" 5.18х1072 1.83x10 
Pr 241x10? 1.77х10% 6.74x10? 2.48x10? 1.03х107 7.58x107 4.11x10? 1.07x10° 1.89x10™ 1.63х1072 5.95x10® 2.11x10 ^! 
Nd 8.86х10° 6.52x109 2.48х10* 9.12x10? 3.80x10? 2.79х10% 1.52х10* 3.95x10? 6.94x10"! 6.01x10? 2.19x10? 7.76x10" 
Sm  L70x10? 125x105 4.77x102 1.75х10° 7.30х10* 5.36x107 2.91х107 7.59x10^ 1.33x10"" 1.16х1072 4.21x10-^ 1.49x107! 
Eu 3.64х10* 2.68x107 1.02x10? 3.75x10^ 1.56x10^ 1.15x107 6.23x10° 1.62x10* 2.85х10:2 2.47х1077 9.01x10-^ 3.19х107° 
Gd  1.68x10? 123x10$ 4.70x10° 1.73x102 7.19х107 5.29х107 2.87x10? 7.48х107 1.31x10"" 1.14х1072 4.15x10® 1.47x10" 
Tb 2.23х10* L64x107 623x10$ 2.29х107 9.54х10° 7.01x10? 3.81x10 5 9.93x10? 1.74x10 ? 1.51х1077 5.51x10 ^ 1.95х1077 
Dy  L17x10? 8.60x107 3.27х107 1.20x10? 5.01х107 3.69x107 2.00x10? 5.22x10* 9.17х107 7.94x10? 2.89х107% 1.02x10" 
Ho  226x10^ 1.66х107 6.34х10% 2.33х107 970x103 7.13х10* 3.87х10% 1.01х107 1.77x10-? 1.54х1077 5.60x10-4 1.98x10 
Er 7.19х10* 529x107 2.01x10° 7.0х107 3.08х10* 2.27x107 1.23x10? 3.21х107 5.64x10-? 4.88х1077 1.78х1077 6.30x10-? 
Tm  9.08x10? 6.68х10* 2.54х10% 9.34x10? 3.89х107 2.86x10? 1.55х10% 4.05x10? 7.12х1077 6.16x10™ 2.25x10™ 7.96х10:5 
Yb 5.99х10* 440x107 1.68x10? 6.16x10^ 2.57х10* 1.89x107 1.02х10° 2.67х107 4.69x10 ? 4.06х1077:1.48х1077 5.25х1077 
Lu  887x10? 6.52х10* 2.48x10° 9.12x10? 3.80x10? 2.79x10* 1.52х10% 3.95х107 6.95х1077 6.02x10^^ 2.19х1075 7.77х1079 
Y 5.99х10° 440x105 1.68х107 6.16x10? 2.57x10? 1.89x10° 1.02x10* 2.67x10? 4.69х107" 4.06x10'? 1.48х1077 524x10 "! 
Sc  345x10? 2.54х10% 9.67х10° 3.55x10? 1.48x10? 1.09х10% 5.91x10? 1.54x10? 2.71x10"! 2.34х1072 8.55х1077 3.03x10™ 


Total 5.87x10? 4.32x10? 1.64x10? 6.04x10? 2.52x10? 1.85x10? 1.00x10? 2.62x10? 4.60x10"? 3.98x10"! 1.45x10™ 5.14x10 1? 


Note: REE, rare earth element; La, lanthanum; Ce, cerium; Pr, praseodymium; Nd, neodymium; Sm, samarium; Eu, europium; Gd, 


gadolinium; Tb, terbium; Dy, dysprosium; Ho, holmium; Er, erbium; Tm, thulium; Yb, ytterbium; Lu, lutetium; Y, yttrium; Sc, 
scandium. 
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Table S2 Non-carcinogenic and carcinogenic risk for REEs in plants of the northeastern Qinghai-Tibet Plateau, 
China 


Non-carcinogenic risk for children Non-carcinogenic risk for adults Carcinogenic risks 
REE 
Ingestion Inhalation Deine Total Ingestion Inhalation pe Total Ingestion Inhalation Dermal Total 
adsorption adsorption adsorption 


La 5.80x10^ 427x107 1.62x10? 5.97x10'^ 2.49х10* 1.83x10 7 9.92х10% 2.59х10 4.55x10- 3.94x10® 1.44x10® 5.08x10-? 
Ce 124x10? 9.14х10 3.48x10? 1.28x10? 5.33x10* 3.92х107 2.13х105 5.54x10 ^ 9.74x10^? 8.44x10-? 3.08x10? 1.0910" 
Pr 128x10^ 938x105 3.57x10° 1.31х107 5.47x10? 4.02х10°* 3.57х107 5.83х107 1.00x10-? 8.66x10-^ 3.16x107^ 1.12x10-? 
Nd 471x104 346x107 1.32x10? 4.84x10'^ 2.02х107 1.48х107 8.05х10% 2.10x107 3.69х1077 3.19х1077 1.16x10 ° 4.12x10? 
Sm 95x10? 6.80x10* 2.59х10% 9.52х107 3.96x10? 2.91х10* 1.58x10'$ 4.13x10? 7.25x10 ? 6.28x10 ^ 2.29x10 8.11x10" 
Eu 222x10? 1.63x10* 621x107 2.28х107 9.50х107 6.98x10° 3.79x107 9.88x10°° 1.74x10 ? 1.50x10 ^ 5.48x10 ^. 1.94x10. ^ 
Gd 924x105 6.80x10% 2.59х10% 9.51х107 3.96х107 2.91x10* 1.58x10 5 4.12x10? 7.24x10 9 6.27х1075 2.29x10™ 8.10х1077 
Tb 123x103 9.06x10? 3.45x107 1.27х107 5.28x10'5 3.88x10? 2.11х107 5.49x10° 9.65x10-^ 8.36x10-5 3.05x10-5 1.08x10? 
Dy 6.64x10? 4.88x10% 1.86x10° 6.83x10° 2.84х107 2.09x10* 1.13x10 5 2.96x10°% 5.20х10:7 4.50x10^^ 1.64x10™ 5.81х1077 
Ho 131x102 9.65x10? 3.67x107 1.35x10? 5.62x10 5 4.13x10? 224x107 5.85x10 5 1.03х10:7 8.90х1077 3.25x10 5 1.15107 
Er 3.75x10? 2.76x10? 1.05х10% 3.86x10? 1.61x10? 1.18x10* 642x107 1.67x10? 2.94x10 P 2.55x10^ 928x105 3.29х107° 
Tm 526x105 3.86x10? 1.47x107 5.41х10% 2.25х10% 1.66x10° 8.99х107 2.34x10 5 4.12x10 3,57x10 5 1.30x10 5 4.61х107'* 
Yb 348x102 2.56x10? 974x107 3.58x10? 1.49x10 1.10х107 5.95x107 1.55x10°% 2.73x10 9 2.36x10 8.61х1075 3.05x1073 
Lu 527x105 3.88x10? 1.48х107 5.4210 2.26x10'5 1.66x10? 9.01х107 2.35x10° 4.13x10- 3.58x10-5 1.30x1075 4.62х107* 
Y 3.51х10* 2.58x107 9.82х10% 3.61х107 1.50х10'* 1.11х10 6.00x10° 1.56x104 2.75x10-? 2.38x10-? 8.68х107'* 3.07х107° 
Sc 155x10* 1.14x107 4.34х10% 1.59х107 6.64x10? 4.88х107 2.65x10°% 6.91x10? 1.21х1077 1.05х1077 3.83х107'* 1.36х107° 


Total 3.31x10? 2.43х10% 9.26x10? 3.40x10? 1.42x10? 1.04x10% 5.80х107 1.48х10 2.59х10:7 2.25х1077 8.19x10 ^? 2.90x10"!! 


